Electronic states in a magnetic quantum-dot molecule: phase transitions and 

spontaneous symmetry breaking 
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We show that a double quantum-dot system made of diluted magnetic semiconductor behaves 
unlike usual molecules. In a semiconductor double quantum dot or in a diatomic molecule, the 
ground state of a single carrier is described by a symmetric orbital. In a magnetic material molecule, 
new ground states with broken symmetry can appear due the competition between the tunnelling and 
magnetic polaron energy. With decreasing temperature, the ground state changes from the normal 
symmetric state to a state with spontaneously broken symmetry. Interestingly, the symmetry of a 
magnetic molecule is recovered at very low temperatures. A magnetic double quantum dot with 
broken-symmetry phases can be used a voltage-controlled nanoscale memory cell. 
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In diluted magnetic semiconductors, the spins of static 
impurities interact through mobile carriers and can form 
a ferromagnetic state [lj, |2j . Since the carrier density in 
semiconductors is a voltage-tunable parameter, the fer- 
romagnetic state also becomes controlled by the voltage 
[3- For technology, this may give an important advan- 
tage compared to the convectional memories based on fer- 
romagnetic metals and controlled by the magnetic field. 
When a single confined carrier interacts with magnetic 
impurities, it forms a new stable state, called a magnetic 
polaron |4j. One important class of confined nanostruc- 
tures is the quantum dots (QDs) where the number of 
trapped carries can be easily changed by the voltage ap- 
plied to the top gate ||- Information in a single QD 
made of magnetic semiconductor can be stored in the 
form of spin polarization and therefore such a QD can be 
viewed as a nanoscale magnetic memory element. Mag- 
netic QDs as memory elements have some properties, 
that look very attractive for the technology: (1) small 
sizes, (2) small number of carriers, and (3) voltage con- 
trol of the number of carriers. Currently, the physics of 
electronic states in magnetic QDs is an active field of re- 
search [1, H S S EH ■ Here we make a logical step from a 
single magnetic nanocrystals toward QD molecules and 
show that a magnetic QD pair has unique physical prop- 
erties that may also be useful for device applications. 

In this letter, we consider a magnetic double QD with 
one hole. Using the mean field theory, we calculate the 
physical properties of magnetic polarons formed due to 
the Mn-hole exchange interaction. We find that this sys- 
tem undergoes two phase transitions (Fig. [IJ. At high 
temperatures, the kinetic tunnelling energy of the hole is 
larger than the magnetic polaron binding energy and a 
symmetric state (with equal hole probabilities for both 
dots) is realized. With decreasing temperature (T), the 
local magnetic energy becomes large enough to trap the 
hole in one of the dots. This self-trapped process sponta- 
neously breaks the symmetry of the system. At very low 



temperature, the local Mn-spin polarization can become 
so strong that the symmetry is recovered. Our broken- 
symmetry polaron has several unique properties: it ap- 
pears in the confined geometry, is voltage-controllable, 
and may vanish at very low T. The above phase transi- 
tions occur also with changing the energy barrier between 
the dots. Since a double QD can be controlled by the gate 
voltage, the above phases can be prepared and read using 
electrical means. 

Model. The Hamiltonian of the system composed of 
two magnetic QDs and one hole has a form 



hh 



p2 



U (R) 



ti -S 

^Jz Jzj 



(1) 



where R = (r, z), r and z are the in-plane vector and the 
vertical coordinate, respectively; j z is the z-component 
of the hole spin and j z = ±3/2; the operator S z = 
Si Z S(R — Ri) is the spin operator associated with the 
Mn subsystem; here Si^ z and Ri are the single-impurity 
spin and position, respectively; i is the Mn-impurity in- 
dex. In this study, we will model the in-plane hole motion 
by a parabolic potential and the vertical motion with a 
square double well i.e. ?7o(R) = u i z ) + TO h/i w o?" 2 /2) 
where u{z) is a double well potential in the z-direction 
(Fig. and ujq is the in-plane frequency, rridh — m * is 
the mass of heavy holes. In our model describing disk- 
like self-assembled QDs, the light-hole states are assumed 
to be split from the lowest heavy-hole states. In the ab- 
sence of Mn impurities, the hole wavefunction can be 
written as tpj = (j}j(z)(p n!C ,n y (r), where J = (j,n x ,n y ), 
n x(y) — 0,1,2..., and j — (s,a,k z ) label the discrete (s 
and a) and continuous states (k z ) related to the z-motion. 
Above we assumed that our QD pair has only two bound 
states, symmetric (s) and antisymmetric (a). 

Within the mean field theory, the Hamiltonian be- 
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FIG. 1: (a) The potential u(z). (b) Schematics of the mag- 
netic QD molecule and three magnetic-polaron phases with a 
single hole. 



comes 



H, 



hh 



2m* 



+ U (R) - QxMnNojz ■ S z , 



(2) 



where XMn is the reduced Mn concentration, Nq is the 
number of cations per unit volume, and «5 Z (R) is the local 
average Mn spin 



(3) 



where S=5/2, B$ is the Brillouin function, To describes 
the effect of anti-ferromagnetic interaction between Mn 
impurities, j z (R) =< ^(B/)|j,*(R - R')|V>(R') >, 
where -0 is the hole wave function. 

In the spirit of mean field theory, the hole moves in the 
presence of the effective spin-dependent potential, 



U = U + SU 7 SU(R) = -$x Mn (R)No3, ■ S 



(4) 



The corresponding ground state wave function has the 
form: rp — ipo\ T>, where | t> is the hole state with 
j z = +3/2. The spatial wave function ipQ should be 
now determined. To the lowest order perturbation theory 
(PT) |l3j , the ground-state wave function can be written 
as = <pG{z)fofi(r), with 



4>g = ac/)s + VT - a 2 cf) a , 



(5) 



where <j> s and 4> a are the bound states in the double 
well. Then, the real parameter a should be deter- 
mined by the variational method. From eqs. [21 El and 



j z (R) = ^ip^,(R) we can obtain the polaron binding en- 
ergy and the total energy 



E b (T) = -%x Mn N S J ^RmfBsi ^^ )], 
El ot {a) = a 2 E s + (1 - a 2 )E a + E b . 

(6) 

The general variational wave function can be written as 
ipG = <fc(z)vo,o( r ) + J2j a J' l Pjj where the second term 
will be considered as perturbation. The zero-order de- 
generate PT provides us with two orthogonal wavefunc- 
tions tpQ and ip% 0]. The corresponding first-order PT 
contribution to the energy is given by eq. Then, the 
second-order correction to the energy becomes 
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where V = {'y,n x ,n v }, 7 = (G,E,k z ), n. 
{0,1,2,...}, T ^ {G, 0,0}, {£,0,0}; here £ G and ^% de- 
scribe the ground state calculated within the lowest-order 
PT and k z labels the delocalized states for the z-motion. 

Phase transition as a function of temperature. 
To find the variational ground state, we calculate the 
total energy E tot (a) = E] ot {a) as a function of a at 
various temperatures (Fig. [2]). The quantity a is the 
hole amplitude related to the symmetric state. The ma- 
terial parameters arc the following: /3Nq = —1.8 eV, 
Nq = 23 nm~ 3 ,Zn = y/ h/m*uo = 2.5 nm, 6 = 4 nm, 
w = 3 nm, Uq — 90 meV, m = 0.38 mo, XMn = 0.005, 
and Rmu = 15.0 nm. This parameter set represents 
a InAs/GaAs QD doped with Mn impurities. One can 
clearly see three different phases. For high T, the ground 
state corresponds to a 2 = 1, i.e., the ground state is 
symmetric. With decreasing T, the ground changes from 
a 2 = 1 to almost a 2 — 0.5. It is a state of hole staying 
mostly in one dot. Thus, the symmetry of the system 
becomes broken. For very low T, the symmetry is recov- 
ered. 

The probability a 2 (Fig. |3J) as a function of the ef- 
fective temperature T* = T + T clearly shows the ex- 
istence of three phases. We may consider the differ- 
ence of hole probabilities of being inside the dots S = 

\Pupper Blower 

I = 2a v / (l-a 2 ) as an order parameter 
(as shown in the inset of Fig. EJ). The order parameter 
5 is nonzero for the broken-symmetry phase at interme- 
diate T* and is zero for the symmetric state at high and 
low T*. 

Basically, the physical picture is the following: the in- 
teraction between the hole and Mn impurities leads to 
the ferromagnetic ordering of Mn-spins inside the double 
QD. The resulting polaron state tends to be as local- 
ized as possible because the magnetic binding energy in- 
creases with decreasing the localization strength of hole. 
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FIG. 2: The total energy of magnetic polaron as a function 
of the wavefunction parameter a at different T* . The arrows 
show the ground states of the system. 

The latter can be seen from eq. the Brillouin function 
strongly increases with increasing the hole wavefunction 
amplitude, ip 2 - Therefore, in our system there is a com- 
petition between the kinetic energy (tunnel splitting) and 
the magnetic energy. At high T, due to the thermal fluc- 
tuations, the interaction between the hole and Mn impu- 
rities is weak and the kinetic energy wins. In this case, 
a symmetric ground state is realized. With decreasing 
T, the Mn-hole interaction becomes stronger and, below 
some critical temperature T C 2, the magnetic localization 
effect overcomes the tunnelling. The hole is now trapped 
in one of the dots since such a spatial configuration low- 
ers the total energy. However, for very low T*, the hole 
is able to polarize the Mn spins inside the QD pair much 
stronger and the magnetic potential SU becomes almost 
flat. Therefore, the symmetric orbital state is recovered. 
There is an obvious condition for the existence of the low- 
T phase: T* x > T Q (T c i > 0). In other words, for the ex- 
istence of the low-T phase, the Mn-hole interaction must 
be strong enough to overcome the antiferromagnetic in- 
teraction. Fig. shows schematically the three phases. 

We should point out that the symmetric polaron at 
low T has a much stronger Mn spin polarization than 
the symmetric state at high T. The temperature T* 2 
corresponds to a smooth, type-II phase transition, while 
the low-T transition is sharp and can be classified as a 
type-I phase transition. The low-T transition occurs as a 
jump between the broken-symmetry phase (a w ±l/v2) 
and the symmetric state (\a\ = 1). 

The temperature T* 2 can be found analytically by ex- 
panding the total energy in the vicinity of the type-II 
phase transition: 



k B T; 2 = P ',°-F^ Sd 3 R W3## ~ (8) 
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FIG. 3: The wavefunction factor a as a function of the effec- 
tive temperature; the function a 2 (T*) shows the appearance 
of broken-symmetry phase (a 2 < 1). Inset: The effective or- 
der parameter (degree of asymmetry) as a function of T* . 
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FIG. 4: Temperatures T* 2 and T*i as functions of the Mn 
concentration. 

Fig.^shows numerically calculated T* 2 (xmu) (dots) and 
the data obtained from eq. |S] (line). Our numerical and 
analytical results agree well. The temperature T* 2 oc 
XMn- This linear dependence may be understood qual- 
itatively. For a higher Mn-consecration, stronger ther- 
mal fluctuations (higher T) are needed to destroy the 
broken-symmetry polaron phase. In contrast to T* 2 , the 
temperature T* x decreases with XMn (approximately as 
l/a^Mn)(Fig- EJ- Qualitatively, such a dependence can 
be explained as follows. The low-T symmetric phase ap- 
pears when the Mn spins are strongly polarized. There- 
fore, a lower T*^ is needed to polarize a larger number of 
Mn spins. 

As discussed above, the transition temperature T* 2 de- 
pends on the competition between the kinetic and mag- 
netic energies and qualitatively corresponds to the condi- 
tion \E a — E s \ ~ SU. T* 2 depends on a material. In Fig.0 
we show the results for a CdMnTe QD. For the parame- 
ters, we use: (3Nq = —1.0 eV, Nq = 15 nm -3 , lo = 4 nm, 
6 = 8 nm, w — 5 nm, Uq = 150 meV, m = 0.1 mo, 
XMn = 0.005, and Rmu — 15.0 nm. The effective tran- 
sition temperature is 14 K and larger than the typical 
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antiferromagnetic CdMnTe temperature To = 3.6 K. 

One may notice that the transition temperature T*^ 
can be very small (Fig. inset). It may be impossi- 
ble to observe the third phase (T* t < To), if some anti- 
ferromagnetic interactions are present. However, due to 
to the tunability of QD systems, we may design a sys- 
tem with a relatively high T* x , Fig. [3] shows the phase 
transitions in a designed system with T* x ~ 2 K . For the 
parameters, we took: (3N = —2.2 eV, N = 10 nm -3 , 
lo = 4 nm, B — 5.3 nm, w = 5 nm, U$ = IbOmeV, 
m — 0.1 mo, iMn = 0.005, and Rmu = 3.0nm. 

In recent experiments [l2|, self-assembled double QDs 
were built into transistor structures with top and back 
contacts. By using applied voltage in these structures 
one can change the number of trapped carriers and also 
modify their wave functions. Two stable magnetic states 
in our system (Fig. QJj) can be considered as a bit of clas- 
sical information. In these states, the hole mostly resides 
either in the upper dot or in the lower dot. Using the gate 
voltage, our system can be prepared in the desired state. 
The following electrical readout of the polaron state can 
be done, for example, using capacitance spectroscopy 
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FIG. 5: Left: Phase transitions in a specially designed mag- 
netic QD pair with a larger T c \. Right: The ground-state 
wave function for different barrier widths for the Mn-doped 
InAs/GaAs QD pair. Inset: phase transition as a function of 
the barrier width. 

Quantum phase transition. We may look at the 
physics from another point of view. The kinetic energy 
of the hole, i.e. the hopping between the dots, is related 
to the width of the barrier between QDs and thus can 
be tuned by appropriate growth process. In Fig. \5\ we 
show the ground-state wavefunctions for various widths 
of the barrier. With increasing of the width, the ground 
states become more asymmetric. It is clear that there is 
a quantum phase transition with the increasing of barrier 
width. The critical width is around 2 nm (see inset of 
Fig. 05 . 

Second order corrections to energy. All the previ- 
ous discussions are based on the perturbation theory ap- 
plied to the nonlinear Schrodinger equation. The second- 
order correction to the energy can be calculated by eq. 
For the sets of parameters used above, it is of a few per- 



cent of the first-order result. We also found numerically 
that, in the vicinity of transition temperatures T* x and 
T* 2 , the second-order corrections are even smaller [l4j . 
This allows us to give the results for T* x and T* 2 in a 
wider range of XMn (see Fig.0J. 

In conclusion, we have studied a magnetic QD molecule 
with a single hole. With decreasing temperature, the 
magnetic polaron ground state undergoes two phase tran- 
sitions. At a higher temperature, the normal symmetric 
state turns into the polaron state with broken symmetry. 
At lower temperatures, a symmetric polaron state can 
be recovered. Our results suggest that a magnetic QD 
molecule can be used as a nanoscale magnetic memory 
cell controlled electrically. 
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